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A statistical theory of the behavior of rays scattered at successive ionospheric reflections 
in propagating around a spherical earth is given. The theory predicts the properties of the 
antipodal area, the deviation of bearing to be expected at places off the main transmitting 
beam and the spread of the incoming fan of rays at all distances. Special properties appear 
to be associated with reception at an angular distance of about 255° from the transmitter. 
The effects of an absorbing region on the received direction and spread of rays are discussed. 

1. Introduction 

In two previous papers [Whale, 1959; Whale, 1962] the effects of scattering by the iono- 
sphere in long-distance shortwave propagation have been discussed. In the former paper the 
concepts of small-circle propagation and the antipodal area were introduced; in the latter a 
statistical treatment of the scattering problem for a flat earth was presented. In the present 
paper, this treatment is modified to cover the case of a spherical earth and thus provides 
solutions to many of the problems which were treated less satisfactorily previously. In par- 
ticular, relations are derived for the mean hearing deviation to be expected when signals ate 
received from a beamed transmitting antenna aimed in some other direction. 

Many of the effects predicted by this theory have been observed experimentally. These 
experimental results will be described in detail in a further paper. 

The assumption is made that the ionosphere is a spherical reflecting and scattering surface 
concentric with the surface of the earth. In practice, this is not quite true; the eccentricity 
of the effective reflecting surface of the ionosphere with respect to the earth's surface leads 
to two effects. The first of these is a shift of the position of the antipodal area away from t he 
point geometrically antipodal to the transmitter. The second, and less obvious effect relates 
to bearing deviations at places about 255° angular distance from the transmitter. At this 
distance, the theory predicts that, under conditions of spherical symmetry and uniform trans- 
missibility of the ionosphere, the mean signal direction should be along the great circle irre- 
spective of the sharpness or direction of the transmitter beam or of the value of the ionospheric 
scattering factor. 

In the last section, the effects of an absorbing region on the received signal are discussed. 

2. Omnidirectional Transmitting Antenna 

Consider an omnidirectional transmitting antenna located at T in figure 1, where the 
great circle to the receiver is T—R. Let TU, UV, VW be successive traces of the ray 
path on the surface of the earth so that 777 is from the ground to the ionosphere, UV is from 
the ionosphere back to the ground and so on. We will assume that the overall effect of scatter- 
ing at the ionosphere and at the ground can be represented by an average effect of the same 
magnitude at all reflections, whether at the ground or the ionosphere. Let TU= UV=VW= 
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Figure 1. Scattering at successive reflection lines 
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. . . =\, where all distances are measured as angles subtended at the earth's center. Let 
the scattering at each reflection be described by the scattering function 



P(0)oc exp-A<t> 2 , 



(1) 



i.e., an incident plane wave in the direction <£=0 becomes a fan of waves, the power in the 
directions contained between the angles <t> and 0+<i<£ being given by (1). The spread of this 
fan of waves can be specified either by the standard deviation, a, where 

1 0.707 



■y/2A -yJ'A 

or by the angle between the direction of the maximum (</>=0) and the half-power point given by 

, 0.832 



Thus 



V3 

a'= 1.18(7, or the beam width = 2a' =2.35<r. 



We can now proceed as follows. Let the primed functions refer to the rays incident at 
one of the scattering reflections and the unprimed functions refer to the emergent rays. The 
power distribution along the line A x is given by 



and 

since 

where 
Then, 



P(Ai)ac exp—A<fi 

AS! 2 
P'(A 2 )oc exp-^S: 
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S 2 = sin A 2 , S=sin X. 
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Performing the same integration on P'(A 3 ) as in (4), we have 



P(A 3 )oc exp— A Q2 , Q l , Q2 0§. (7) 

The general forms are thus: 



s 2 +st+sr 



P'(A w )cc exp-^ g2+/S;+ gff > >+g 2 i «• (8) 

C2 

P(A„)cc exp-^i ^2 +aS 2 +>S 2^_ q^r 2 *J. (9) 

The denominators are relatively slowly varying functions so that we may use the approximation: 

S n — sinA„~sin n\. 

In any actual case, the angular distance from the transmitter to the receiver is known 
(TR=A n ) and X must be deduced from the number of hops (either assumed or found from 
measurements of vertical angles). 

The series in the denominators of (8) and (9) can be summed to yield: 

n , , K . 2A sin A -sin 2 A w ., /1AN 

P'(A„)cc exp — — . — = — — & (10) 

(n— 1) sin X— cos A„ sin (A„— X) 

d/ k \ 2-4 sin X sin 2 A„ a , . 

P(A„)oc exp r-— — : — — 2 - (11) 

/?- sin X— cos (A„,+X) -sin A n 

In a previous publication, an approximate value for the spread of a wave received at a 
distant point was obtained using a theory based on small circle propagation [Whale, 1959]. It 
is interesting to compare this former result with the one obtained above. 

The standard deviation of the spread calculated from (10) can be written as: 

1 , Tr %A sin X shr A„ 

<t= ,,= where K= 



yj2K (n—l) sin X— cos A n sin (A„ — X) 

For comparison, put X=9°, and A such that a=2° at A w =90°, i.e., /1-2000, in this case. 
The small circle formula becomes 

(sin) A0=sin cr-0.035 tan ~ 

z 

The two curves are plotted in figure 2. There are significant differences between the two. 
The statistical theory given above is useful up to any distance, even beyond the antipodal point, 
while the previous theory could not be extended to that distance. 

3. Directional Transmitting Antenna 

The directivity of the transmitting antenna may be introduced as a term involving 6 
(fig. 1). If this directivity is of the form 

P(d)ocexv-Bd 2 , (12) 

then 

P(Ai)oc exp-Bdl .exp-Atf (13) 

where the point at which the observations are made is specified in terms of the coordinates 
(A 2 , $i) i.e., colatitude and longitude measured from the center of the antenna beam. 

If we proceed as in section 2, the expressions for the distribution of directions at greater 
distances rapidly become very complicated. The general expressions are given below. If the 
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distributions along successive lines are required, it is more economical of time to compute each 
distribution from the previous one using the recurrence relations. If the expression forP'(A„) 
is given by 

P'(A w )oc exp-(L'd 2 n -M'e t £ n +N'ti), (14) 

then 

P(A„)ocexp~|[/y 

If this is now written 
then 



1 2 AM' ajL AN' ±2 \ ( . 

V)\ 6l ~A+N' 6 ^ + A+N' *r (15) 



(16) 



4(A+N* 

P(A Joe exp-iLdl-Mdn^+Nc^l), 
P'(A n+1 )oc exp- ^Lel +1 -^2L £+M ^] *„+&.+, 

[ i (£)*^£X§fM% ! )>«} m 



+ 



where 



<S'=sin X. S„=sm nX. 



The general forms of the expressions for the probability distributions may be written as 
follows: 



If P(A„) oc exp-(L9l-M9 n <t> n +N<t>l) 
then 



(18) 
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Similarly, if /"(A„) oc exp-CL'^-MXln+iV'^), (22) 

then 
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(25) 



The series appearing in (19) to (25) maybe summed but the relations between the various 
quantities can be seen rather more readily when they are written in the above forms. 

At any place specified by the coordinates (0 n , A n ) the spread of the incoming wave may be 
found directly from the expression for P'(A B ). On the antenna axis (i.e., in the direction of 
the main beam) 6 n =0, so that 

P'(A n )«exp-iV'£ n 2 . 



The spread of the directions is thus given by 

1 



\^V 



(26) 



At places off the axis of the main beam, the received signals are no longer symmetrically 
spread about the great circle direction to the transmitter. The direction from which maximum 
signals are received is given by dP'/d£=Q i.e., 

Z=2W' (20 

To obtain the polar diagram of the transmitting antenna as measured with an omni- 
directional antenna at the receiving site, the total power received as 6 is varied is found from 

r P'^i.e., 

c/ — 00 

PWcxexp-(/y-^)^ (23) 

4. Spread of Incoming Directions 

The directions of the signals arriving at any place with coordinates (A, 0) are contained in 
the expression 

P , (A)oc exp-(Z'0 2 -M'0£+A^ 2 ). (29) 

On the antenna axis 6=0, so that cr= , • This quantity is plotted in figure 3 using 

the values of the constants X=ll°, .4=1000, B=0, 10, 30, 1000. 
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The curves in figure 3 may be used for values of A other than 1000. If the desired value of 
A = a, then the spread read from the curves must be multiplied by ^{lOOO/a} while B is new 
the ratio 10005/a, 

The same spread is obtained for places the same distance from the transmitting antenna 
but off the main axis of the beam. The maximum of the incoming directions is, however, then 
no longer along the great circle connecting the transmitter to the receiver (except in one special 
case treated later on). 

A noticeable feature of the curves in figure 3 is that, at an angular distance of about 255° 
from the transmitter, the spread of the incoming directions is apparently independent of the 
beam width of the transmitting antenna and thus depends only on the reflection scattering 
coefficient. 

The value of the standard deviation of the spread at this particular distance is given by 



1 ln—1 

TVT 



where n=number of steps. 



(30) 



For X=ll°, n=23 at 253°, so that <t^4.°2 (A=1000). 

The vicinity of A=255° has further interesting properties, as will be seen in the next 
section. 

At the antipodal point (A=180°) the value of N' can be found, by summing the series in 

AB 
(25), to be given by N'= — r— 

A+^B 



Since the standard deviation of the spread of power in the transmitter beam is given by 

1 
y2£ 

the ratio of the spreads of the received fan of waves at the antipodal point to the transmitted 
fan is given by 
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Figure 3. Spread of incoming directions as a func- 
tion of the distance from the transmitter for vari- 
ous transmitting antenna beamwidths. 

The half-power beam-width of the transmitting antenna is zh—n deg. 
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5. Deviation of Incoming Directions From the Great Circle Direction 

The maximum of (14) is obtained at 



sin six 



M' 

~~2N' 



If the receiver is on the main axis of the transmitter beam, 0=0, so that the received direc- 
tions have a maximum along the great circle to the transmitter. Off the main axis, the ratio 
of the deviation of the maximum of the received directions to the angle between the trans- 
mitter beaming direction and the great circle from the transmitter to (lie receiver is given by 



in 



M' 



e T 2N' 



(31) 



This quantity is plotted in figure 4. 

At the antipodal point to the transmitter, the received angle follows the transmitted angle 
exactly, as would be expected. There is a further interesting region, distant about 257° from 
the transmitter, where the received bearings are always clustered about the great circle direc- 
tion irrespective of the bearing direction and beam width of the transmitter and independent 
of the scattering coefficient A. From the expressions for M' and N f in (24) and (25) it can be 
seen that (31) goes to zero when the bracketed term in (24) goes to zero. This term docs not 
contain A or B. 

The line at about 255° from the transmitter where the spread is independent of the trans- 
mitting antenna beam width and where the mean direction of the incoming waves is always 
along the great circle direction, could perhaps be called the "reciprocal line" to the trans- 
mitter. These conditions hold for an ionosphere which is concentric with the earth and in 
which all paths are possible. It is obvious that experimental measurements taken at this 
distance offer a unique opportunity to differentiate between scattering effects and other iono- 
spheric effects. 

The sign of the variation of £ max with T is important. From figure 4 !, if the receiver is at 
W with the transmitter beam aimed along the line TH, then 6, the 4 longitude of W with respect 
to TR is positive and £ 3 , the deviation, is positive when measured in the direction shown. 
At receiving points more than 180° distant from the transmitter, £ ma , must be measured in a 
clockwise direction from the great circle (long path) connecting the receiver to the transmitter. 
The reciprocal line effect discussed above may be difficult to observe since the short path to 
the transmitter is only about 100° of angular distance long. Backward radiation will then, of 
course, lead to the bearing deviations expected at Fin figure 4. 



Figure 4. Ratio of the deviation of the direction of 
the maximum of the incoming fan of rays from the 
great circle direction to the angle between the trans- 
mitter beam and the great circle direction to the 
receiver. 
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6. Effects of an Absorbing Region 

The presence of a region in the ionosphere which inhibits the propagation of radio waves 
(for example, the absorbing regions near the auroral zones or a region of low critical frequency) 
may be introduced into the analysis by a suitable modification of the distribution along any 
line A. 

Let the undisturbed distribution of signals emergent from the line A n be given, as before, 
by the expression 

P(A n )oceW-(Ld 2 n -Md n <t> n +N<t> 2 n ). (32) 

Consider an absorbing region placed as shown at E in figure 5a. This is along one of the 
latitude lines in the coordinate system referred to T as a pole. Then the shaded areas are the 
shadow zones as far as great circle propagation is concerned; APis the point antipodal to the 
transmitter. The edge E is at 0=0, i.e., at the center of the main beam of the transmitting 
antenna. 
Then 

o 
P'(A„+i) oc exp— {L'6\ +l — M'e„ +1 t n+l +N'& +1 ) with £ now limited to the values £>-^ 8„ +i . 



(33) 



Thus 

P(A, i+I ) oc J Sn ex V -A(<j>-£) 2 -P' (A„ +1 ><£ oc exp-(Z0 2 n+1 -M0 B+ A + , 

+ « t ,,[i + e rf {vzw(¥-^)}) 

Since we are mainly interested in regions well within the shadow of the stop, we may 
replace the error function by its approximation 

erf( „ )s=1 _ e _5Pl^l 2 . 

v 
This gives, for v greater than zero, i.e., within the shadow area, 

P(A n+l ) ocexp- (Ld 2 -Mdct>+N4> 2 ) -[l + erf(v)]~2 ex^-(Ld 2 -Md<t>+N<t> 2 ). (34) 




Figure 5. Shadow zones created by an opaque stop. 
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The exponential term represents a fan of rays centered on the edge of the stop, with the 

same spread as the incoming beam. The situation can be represented physically as in figure 6. 

The distribution emergent from the line A n is given by (32), and the distribution incident 

S 
on A w+1 is given by (33), provided that £>-<t0. For any 0, the maximum value of the incident 



distribution is given when %=-£$, i.e., from (33), 

o 

P'(A„ +1 ) oc e X p-|i^|'M^+(^)V | 2 

=exp-N (^p) 2 2 > where N is the factor in P(A„) (35) 

For any 0, P'(A W+1 ) given by (33) yields the same spread lor £ so that typical plots for 
P'(A„ +1 ) at different positions within the shadow zone are as shown in figure 6a. This spread 
of £, since it is independent of d, is obtained by putting 0=0 in (33), i.e., P'(A W+1 ) oc exp — N'£ 2 

= exp— N'( -£ yd 2 when £ is written in terms of 6. 
From (35), P' (A w+1 ) aexp-A^/^V, for the value of P'(f) max 

=exp-iV w (%^) 2 * 2 , in terms of f . 

On emerging from the scattering line A n+i , the distributions shown in figure 6a are broad- 
ened by the scattering factor A. Typical shapes would then be as in figure 6b. This is in 
agreement with the expression for P(A n+1 ) given in (34). Thus, for receiving points within 
the shadow of the stop, we may replace the transmitter by a new source at the (n\gc of the stop 
with its polar diagram alined in the direction of the great circle connecting the transmitter to 
the stop edge. The required virtual source at E has a polar diagram of the form 



P(4>)ocexp-7v> 2 . 



(36) 



The analysis given so far in this section certainly applies to receiving points within the 
first shadow zone of the stop, i.e., places like if in figures 5a and 5b. The situation at places 
like J beyond the antipodal point to the transmitter (AP) may be slightly different. Pro- 
vided the distance EJ is less than 180°, the angle which the great circle through E and J makes 



p'(A n+l ) 



Figure 6. Distributions of directions near the edge 
of an opaque stop. 
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exp-N(^±i#) 




with the line T—E—AP is such that the required initial scattering direction is away from the 
immediate shadow zone. In this case the beam-width of the source at the shadow edge is 
given by 

P(4>) oc exp-A<t> 2 . (37) 

These correspond to the very narrow beams shown in the left-hand parts of the distribu- 
tions sketched in figure 6. 

This result may be obtained directly from (33) since, under these conditions we obtain 

P(A n +i) ozexp-(Ld 2 n+l -M6 n+l ct> n+1 +N^ +l )- [1-erf (v)] 



oc 



r~ -J =p- - exp- f(u+A+N' %) -M f -A d 2 -2A % 6<t>+A^\ 

\2(A+N')¥-M'\e-2M lV V ' J 

L ^ J ( 38 ) 

oc exp— A<f> 2 , for large </>. 

7. Conclusion 

The theory presented above is in agreement with many of the observational results which 
have been published previously [Whale, 1956; Whale, 1959; Whale, 1962]. It leads to a much 
better appreciation of the characteristics of the antipodal area than has been possible hereto- 
fore and suggests interesting phenomena occurring at a distance of about 255° from the trans- 
mitter. 

Experimental investigations of the phenomena described above are being conducted; the 
results will be submitted to this journal. 
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